Photon emission from out of equilibrium dissipative parton plasma 
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Using the second order Israel-Stewart hydrodynamics we discuss the effect of viscosity on photon 
production in a parton plasma created in relativistic heavy ion collisions. We find that photon pro- 
duction rates can enhance by several factors due to the viscous effect in a chemically nonequilibrated 
plasma. 

PACS numbers: 



A strongly coupled quark-gluon plasma (sGGP) or a matter in a perfect fluid state is widely expected to be produced 
in recent Relativistic Heavy Ion Collider (RHIC) experiments. Recent measurements has shown that the matter flows 
very rapidly at the time of its breakup into the freely streamiiig hadronic matter. Also the measurements of the 
elliptical flow parameter V2 show a strong collectivity in the flow [J, 0, This would imply that the QGP can have a 
very low shear viscous stress and the ratio of its shear viscosity rj to the entropy density s i.e. ?//s should not be much 
larger than the lower bound l/47r Q. This led to a conjecture that the QGP formed at RHIC is the most perfect-fluid 
found in nature There has been a lot of attempts to determine the viscosity of sQGP 0, d, 0, H, H, The 
first order theory of viscous hydrodynamics is known to give unphysical results. For example, when the Navier-Stokes 
equations were applied to the one dimensional boost invariant expanding flow [11], one finds the expression for the 
temperature to be 
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This describes reheating of the flow and T has a maximum at time r^ax = '''o \ \ + f '^°2° ) ■ This is an un- 
physical behaviour and the first order viscous dynamics is known to have such problems [l^. Il3|. The second order 
hydrodynamics approach developed in the spirit of Israel and Stewarts |1J] removes such artifacts. The second or- 
der viscous hydrodynamics later developed and applied in the context of heavy-ions collisions by several authors 

[H d, [O, [il, [H, iH, HI] . 

It would be interesting to study the role that viscosity can play on the plasma signals. Hard photons are one 
such promising source that can provide information about the thermodynamical state of the plasma at time of their 
production. The plasma created in the heavy-ion collisions is expected to be in a state of chemical nonequilibrium. 
The photon emission from such a plasma has been studied within the framework of ideal hydrodynamics by earlier 
workers [l^, US, HB, UHl ■ In this paper we study the photon production using causal hydrodynamics of Israel-Stewart 



In the center of the fireball in a nuclear collision the viscous stress-energy tensor in the local comoving frame has 
the form [11, [IS, El: 
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with the transverse and longitudinal pressure 



(2) 



Here P denotes the (isotropic) pressure in thermal equilibrium, $ denotes the non-equilibrium contributions to the 
pressure coming from shear stress. We ignore the bulk viscosity in the relativistic limit when the equation of state 
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p = e/3 is obeyed [2^. However, the bulk viscosity can be important near the critical temperature [30, IMI- The 
shear tensor in that frame takes the form tt*-' — diag($/2, <l'/2, — $) consistent with the symmetries in the transverse 
directions. 

To describe evolution of the energy density and the viscous stress (jj we use second order dissipative hydrodynamics 
of Israel-Stewart [H [H [H, [H : 

de 1 

^ = --(e + P-*), (3) 
^ = -^-U('- + ^T^(^-1)) + '-^ (4) 

where f32 = 9/(4e) and Tt^ = 2/?2?7 denotes the relaxation time. Equations (3-4) are written in the local rest frame 
using hydrodynamic velocity u'^ = i(i,0, 0,z), where r = \Jt'^ — W^- Equation of state is required to solve these 
equations. We use ultra-relativistic equation of state : P = \t. 

To describe the chemical non-equilibration while maintaining the kinetic equilibrium, one can use the parton 
distribution 13611 . 



where, u'^ is the four-velocity of the local comoving reference frame. The temperature T is a time-dependent quantity 
and the distribution is multiplied by time and another dependent quantity called fugacity Xq^g (r) to describe deviations 
from the chemical equilibrium. The fugacity parameter become unity when the chemical-equilibrium is reached and 
in general it has the range < Xq^g < 1. The scattering processes gg <-> ggg and gg ^ qq give the most dominant 
mechanism for the chemical equilibration. The master equations describing evolution the parton density are given by 

2 ~ ? \ 




d.inqun = -a^nill-^ , (7) 



where ni{i = q,g) is parton density with unit fugacity 37] and (72 = {'^{gg ^ Qq)) £ind a3{a{gg ^ ggg)) are thermally 
averaged scattering cross sections. It should be noted here that when equation for Ug and nq are added one gets the 

total number density n and the term with \cr2'n'^g ( 1 — jpj^ ) will drop out. This is because due to the the scattering 

process gg qq loss in the gluon density is equal to the gain in quark density and vice verse, 
e and n can be calculated using equation (5) as given below 

n= (Agai +A,6i)r3, e ^ {Xga2 + Xqb2)T^ (8) 

where oi = 16f(3)/7r2, a2 = Stt^/IS for the gluons and bi 9^{3)Nj: /n'^ , 62 = 7T:'^Nf/20 for the quarks. Using 
equations (3-4,6-8) following evolution equations for T, Xq g and $ can be obtained 
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where, the rates R2 = 0.24:NfalXgTln{5.5/Xg) and R3 = 2.1a2T(2Ag - A^)!/^ are defined as in Ref. [H, We 
would like to note that our gluon fugacity equation pTI) differs from that given in Ref. [H, [s^l by a fetor of two 
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in second term in right hand side. We beheve this is a typographical error. In equation Q the first term on left 
hand side is due to expansion of the plasma, while on the right hand side the first term describes effect of chemical 
nonequilibrium and second term is due to the presence of (causal) viscosity. The last term in parenthesis of equation 
PH)) arises because of the chemical nonequilibrium process. It should be noted that equation (9) differ from that 
considered in Ref. [38l|. In their treatment first order viscous hydrodynamics is used which does not require time 
evolution of However such treatment give unphysical results like reheating artifact as mentioned before. 

Elastic (gg <-> gg) as well as nonelastic processes like gg <-> ggg can contribute to the shear viscosity. Shear viscosity 
coefficient was recently calculated for the inelastic process in the presence of chemical nonequilibrium in Ref. (ssj . It 



was shown that r]/T^ ~ Ug/T^ 
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From this one can write [3 
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It ought to be mentioned that this viscosity prescription was not considered considered in Ref. [3^. Kinetic theory 
without invoking nonequilibrium process gives = 3/2t:T. 

Real photons are produced from the annihilation of a quark-antiquark pair into a photon and a gluon (qq gj) 
and by absorption of a gluon by a quark emitting a photon (qg — > 57). Another source of photon production could be 
the bremsstrahlung but its effect can be ignored in the lowest order of a perturbation theory. In order to compute the 
photon production rates one needs to know the underlying amplitude Ai of the basic process involving the annihilation 
or Compton scattering process and the parton distribution functions given by [1^ |40| 
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Here Pi and P2 are the 4-momenta of the incoming partons, P3 of the outgoing parton, and P of the pro- 
duced photon. In equilibrium, the distribution functions ni(Ei) are given by the Bose-Einstein distribution, 
nsiEi) = l/[exp{Ei/T) — 1], for gluons and by the Fermi-Dirac distribution, npiEi) = l/[exp(£'i/T) + 1], for quarks, 
respectively. The factor (|A^p) is the matrix element of the basic process averaged over the initial states and summed 
over the final states. The indicates the sum over the initial parton states. The fugacity factors can enter equation 
(14) when equation (5) is considered 



ni{E)n2[E){l ± ns,{e)) ^ AiniA2n2(l± X^.n^). 



This is can be rewritten as 



AiniA2n2(l ± A3?i3) = AiA2A3nin2(l ± 713) 
+ AiA2(l - A3)nin2 



(15) 



In carrying out the momentum integration it is useful to introduce a parameter fee to distinguish between soft 
and hard momenta of the quark [4l|. For the hard part of the photon rate following [23| we take = 2m^ — 
0.22g^T^ {Xg + Ag/2), where, is the quark-thermal-mass which can be obtained from zero momentum limit of 
quark self-energy in the high temperature limit. The first term on the right hand side of equation (15) can lead to 
the following photon rate [23| using the Boltzmann distribution functions instead of a quantum mechanical ones: 
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Here a and as are the electromagnetic and the strong interaction coupling constants. The second term in equation 
(15) will give, under the Boltzmann approximation, the following contribution to the photon rate: 
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The total photon production rate 2E ^i^i^ can be obtained by adding equations (16-17), is required to be convoluted 
with the space time evolution of the heavy-ion collision. We define [1^ 

where the times after the maximum overlap of the nuclei are ro and and ti and ynuc is the rapidity of the nuclei. Q 
is the transverse cross-section of the nuclei and p_L is the photon momentum in direction perpendicular to the collision 

axis. For a Au nucleus Q ^ 180/to^. The quantity {^E ^J^^^ is Lorentz invariant and it is evaluated in the local 

rest frame in equation (18). The photon energy in this frame, i.c in the frame comoving with the plasma, can be 
written as p±cosh{y — y'). 
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FIG. 1: Temperature, gluon fugacity and quark fugacity for RHIC and LHC. Solid lines indicate the case with the shear 
viscosity, while the dashed lines correspond to the case without viscosity 

In Figure (1), we have shown T,Xg,Xq as function of time. We have solved the equations (9-12) together with the 
initial conditions at n^o from HIJING Monte Carlo model Which are A° = 0.09, A° = 0.02 and To = 0.57GeV 
for RHIC with T,,o = 0.31/m/c and A° = 0.14, X° = 0.03 and To = OMGeV for LHC with t,so = 0.23/m/c. 
Presence of the causal viscosity decreases the fall of temperature due to expansion and the chemical nonequilibrium. 
However if one considers the first order theory, there can be unphysical instability. Fugacity of gluons and quarks 
increase more slowly due to the presence of the viscosity compared to the cases when no viscous effects were included. 
This is because the chemical equilibration is reached here with falling of the temperature. The temperature can 
decrease due to the expansion and chemical nonequilibration. The lowering of T can help in attaining chemical 
equilibrium and which in turn will increase the rate at which the fugacities reach their equilibrium values. Inclusion 
of the viscosity will slowdown the falling rate of the temperature. Consequently the fugacities will take more time to 
reach their equilibrium values. 
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FIG. 2: (Left panel) Photon rate for different rapidities in LHC (ynuc ~ 8.8). (Right panel) Same with the inclusion of viscosity. 

We plot photon spectra by using equation for T,r in solving equations (9-12). The figures (2-3) compare the 
case without viscosity with the case of finite shear viscosity. 

Figure (2) shows the photon spectra emitted at fixed rapidities as a function of transverse momenta p±. The 
photon flux is normalized with the transverse size of the colliding nuclei((5)- For LHC we take: tq = 0.5/to/c, 
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FIG. 3: (Left panel) Photon rate for different rapidities in RHIC (ynuc=6.0). (Right panel) Same with the inclusion of viscosity. 




1.0 1,5 2,0 2,5 3,0 n-"l 

I,( 1,5 2,(1 



FIG. 4: Photon rate for different rapidities in RHIC (ynuc=6.0) and LHC (ynuc=8.8) with kinetic viscosity. 

Ti — 6.25/m/c and ynuc — 8.8. We use equation (fT3|) for tv in solving equations (9-12). The figure compares the 
case without viscosity with the case of finite shear viscosity [35| . 

Figure (3) shows the comparison similar to that of figure (2) but with a set of initial conditions for RHIC: tq ~ 
0.7 fm/c,Ti = Afm/c and ynuc = 6.0. For as = 0.3, shear viscosity to entropy density ratio r//s ~ 0.29. Figures 
(2-3) show that viscous effects enhance the photon flux by a factor (1.5-2). 

Finally, we compare the photon fluxes calculated using equation (jl3p with the fluxes calculated using the kinetic 
viscosity (tjt = 3/2ttT). Figure (4) shows the photon flux calculated using the kinetic viscosity prescription for LHC 
and RHIC. However, we do not find any significant change in the flux for the results obtained using equation ()13|) . 

In conclusions, we have studied the second order dissipative hydrodynamics with chemical nonequilibration. We 
find that the effect of viscosity enhancing the photon flux by a factor ranging between 1.5-2 for the parameter space 
relevant for LHC and RHIC. Our results are in a broad qualitative agreement with the results obtained in Ref . ^38| 
using the first order theory. We also find that the two viscosity prescriptions with inelastic scattering [s^ and the 
one involving elastic collisions only using kinetic theory give similar results for the photon production rate. 
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